Ensemble Equivalence for Counterion Condensation on a Two Dimensional Charged 

Disc 
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Abstract 

We study the counterion condensation on a two dimensional charged disc in the limit of infinite dilution, and compare the 
energy-temperature relation obtained from the canonical free energy and microcanonical entropy. The microcanonical entropy 
is piecewise linear in energy, and is shown to be concave for all energies. As a result, even though the interactions are long- 
ranged, the energy-temperature relation and hence the counterion condensation transition points are identical in both the 
ensembles. 
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Polyelectrolytes are polymers that release counteri- 
ons into a polar solvent, making the polymer backbone 
charged [ll-y]- They are ubiquitous in biological sys- 
tems - examples include DNA F-actin, microtubles 
and tobacco mosaic virus [5]. On increasing the lin- 
ear charge density of the polymer backbone beyond a 
threshold value, counterions start condensing onto the 
polymer @, |7j. The interactions due to the condensed 
counterions lead to complex phase diagrams [3, 0, H, Q ; 
and phenomena such as collapse of an extended poly- 
electrolyte chain [Iol - fl2l |. and ag greg ation of similarly 
charged polyelectrolyte chains pi Il3l H3 |. This coun- 
terion condensation transition (CCT) has been mostly 
studied in the mean field limit, for idealised systems like 
cylinder jUTEEl, planes [TJSl and spheres QJ], due 
to the analytical intractability of real polyelectrolyte sys- 
tems. In recent work, the partition function of a two di- 
mensional system of counterions around a charged disc 
was evaluated in the limit of infinite dilution [201 ] . The 
theoretically obtained critical points of CCT and the de- 
pendence of energy on temperature matched with results 
from Monte Carlo simulations [ill [l6| . 

Electrostatic interactions are long-ranged, i.e., in d- 
dimensions, they decay slower than r~ d at large distances 
r. A possible consequence of long-range interactions and 
non-additivity of energy is the inequivalence of different 
statistical ensembles [21, 22]. This inequivalence may be 
manifested as negative specific heat or magnetic suscep- 
tibility in the microcanonical ensemble, different order 
of phase transitions, and different critical points [2l| - |23| . 
As systems that show CCT are long-ranged, it would be 
of interest to know if CCT is the same in different ensem- 
bles. However, almost all such models are only studied in 
the canonical ensemble. In this paper, we solve the model 
of CCT in Ref. [2fJ of a system of counterions condens- 
ing onto an oppositely charged disc in two dimensions 
in the microcanonical ensemble to address the question 
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of ensemble equivalence of CCT. We first obtain the ex- 
pression for entropy and show that it is piecewise linear 
in energy. This entropy is then shown to be concave ev- 
erywhere in energy and volume. The CCT temperature 
and the energy-temperature relations obtained in the mi- 
crocanonical ensemble are shown to be identical to those 
obtained from the calculation in the canonical ensemble 
|2Qj. Thus CCT, at least in this solvable case, is shown 
to be equivalent in the canonical and microcanonical en- 
sembles, despite the interactions being long-ranged. 

Examples of non-additive long-ranged systems that 
show ensemble inequivalence include spin systems |24| . 
IHJ, self-gravitating systems [26l - [28j . plasma |29| . l30j |. 
two dimensional hydrodynamic systems [3l|, [33] and 
other model systems [33rl35| . In the last decade, there 
have been various attempts to compare the solutions of 
long-ranged systems obtained in the microcanonical and 
canonical ensembles [22]. However, the number of exactly 
solvable models where one can pinpoint the equivalence 
or inequivalence are limited, since the calculation of mi- 
crocanonical entropy is difficult even for systems with 
short-ranged interactions. The solution presented in this 
paper will add to this list of exactly solvable models. 

Consider a uniformly charged disc of charge q and ra- 
dius a. N counterions, each carrying a charge —q', are 
distributed in the annular region between the charged 
disc and a circular boundary of radius R. Overall charge 
neutrality is achieved by choosing q = Nq'. Let be 
the position of counterion i in a coordinate system with 
origin at the center of the disc. The Hamiltonian of the 
system is 
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(1) 



where r. 

x = qq' 



= | r^l and and ry = — r j I • T ne parameter 
= Nq' 2 measures the strength of the electrostatic 
interaction. The permittivity e has been set equal to 
l/(47r). The thermodynamics of th is sy stem was studied 
in the canonical ensemble in Ref. [20]. We now obtain 
the expression for the partition function following closely 
the steps in Ref. [20| . 
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The Hamiltonian can be rewritten in terms of U{ 
ln(r,/a) as 



y i=i 

In [2 cosh(t 



X 

2N 



u i )-2cos%], (2) 



where cos 9ij — n ■ rj/ (rjrj). This Hamiltonian is analyt- 
ically intractable. However, some simplifications occur 
when the limit of infinite dilution, N/R 2 — > 0, is consid- 
ered. This corresponds to the limit of finite number of 
counterions in infinite volume, similar to the limit consid- 
ered in Manning condensation Q. In this limit, counte- 
rions are far away from each other. Thus, \ui — Uj\ 3> 1, 
and the summand in the second term of Eq. @ may 
be approximated by \ui — uj\, since 2cosh(ii; — Uj) sa 
exp — itj|). In this limit, the partition function of 



the system was calculated in Ref. [2C 
slightly different notation), 



and is given by (in 



N N 

z =n n 

fc=0 m = 
m =/= k 



2 (k — m) x 



x exp 



-2(N-k) X L fp_ H(N,k) 



HN,k) V X 
where L = ln(i?/a), /3 is the inverse temperature and 

mm) 



(3) 



2N - k - m + 1 



(4) 



In the limit L — ¥ oo keeping N fixed, the sum in Eq. ^ 
is dominated by the largest summand. Given a value of 
f$Xi the k corresponding to the largest summand is ob- 
tained by solving §x = N/ (N — k + 1) for k and taking 
the integer part [20|. As j3x is increased, k changes by 



unity, resulting in non-analytic behaviour of the free en- 
ergy. These critical points occur at 



N 



N 



= £(fc-l,fc), k= 1,2,... N. (5) 



The first transition is at fix = 1 (k = I), corresponding to 
the Manning condensation temperature on cylinders [6[. 

To show the ensemble equivalence, we now compute 
the entropy in the microcanonical ensemble from the 
canonical partition function in Eq. ([3]). We then show 
that the entropy is concave everywhere and the CCT 
temperatures coincide in both the ensembles. 

The density of states g(E) is defined as 



i f N 

g{E) = w. J Ii driridei 6 ( E ~ H ) 



(6) 



g(E) is obtained from the partition function Z by per- 
forming an inverse Laplace transform, i.e., 



9(E) 



dp Z(/3) exp((3E), 



(7) 



2ni J_ ioo+C 

where the constant c is to be chosen such that the path 
of integration lies to the right of all the poles of Z(j3). 
This corresponds to £(k,m)/x < c for all k and m. We 
evaluate the above integral by the method of residues, by 
closing the contour of integration with a semi-circle that 
is either to the left or right such that the contribution to 
the integral from the semi-circle is zero. If E and k do 
not satisfy the condition, 



E 



2{N-k) X L 
£(k,N) 



>0, 



(8) 



then the contour of integration is closed to the right in 
the complex plane. The closed contour does not enclose 
any of the the poles of the partition function, and, hence, 
the contribution to g(E) from such fc is zero. On the 
other hand, if E and k satisfy the condition in Eq JHJ, 
the contour is closed to the left, and g(E) is the sum of 
the residues of the partition function, and is given by 
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k=k* i = o 

/ ^ k 



2(k-l) X 



exp 



~(N — k)(N — l)£ (k, I) L £(k,l)E 
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m = 
m ^ k 
m ^= I 
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(k-m) (l — m)£ (k, l)_ 



(9) 



where k* is the smallest value of k, given E, that satisfies 
Eq. ([8]). In writing Eq. (j9)), we have used the relation 

£(l,k) -£(m,k) = — , (10) 

derived from Eq. The terms in the summation with 



I > k* do not contribute to g(E) because the summand is 
antisymmetric in k and I. This leads to an upper bound 
k* — 1 for I. Also, the maximum possible value of the 
energy of the system is (N + l)xL, corresponding to ui — 
L for all i. It then follows from Eq. ([5]) that the lowest 
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FIG. 1. The variation of microcanonical entropy S with en- 
ergy E is shown for different L. The data are for N = 3 and 
X = 3. The solid black lines correspond to the limiting curve 
obtained from Eq. <|11[) . They have been shifted upwards and 
extended to the left and right for clarity since the unshifted 
lines are indistinguishable from the curve for L — 200. 



value of k* is 1. 

The density of states g(E) in Eq (JSJ is a sum of expo- 
nentials of the form exp[/(fc, I, E/L)L], where the func- 
tion / is derivable from Eq (|9]). In the limit L — > oo, 
keeping N fixed, the sum is dominated by that k, I which 
maximize the function f.Eis implicitly dependent on L 
and we make this dependence explicit by converting the 
inequality in Eq. (JSJ into an equality by replacing k by 
k* —S, where < 6 < 1. Thus, eliminating E, / becomes 
a function of k, I and k*. It is then easy to show that 
/ increases with I for k > k* , and decreases with k for 
I < k* . As a result, the dominant contribution to g(E) 
comes from the term with k = k* and I = k* — 1. 

By considering only the largest term, we obtain the 
entropy S — \n[g(E)] to be 

lim A = 1 -^l + ^ fc *- 1 )A, 

where f(fc*,fc*-l) = N/(N—k* + l). In Fig. [Q we show 
the dependence of entropy on energy, obtained by con- 
sidering all the terms in g(E) [see Eq. (|S}], and compare 
it with the entropy in Eq. ([TTj) . On increasing L, S/(NL) 
approaches the limiting curve with TV linear portions, de- 
termined by Eq. pT|) . 

We now show that the entropy is concave in energy 
E. In Eq. (TIT]) , k* , though a function of E, takes on 
integer values and is a constant over a range of E. As 
E increases from to its maximum value, k* decreases 
from TV to 1, in steps of unity. Thus, the entropy curve 
consists of N linear segments of slope £(k*,k* — 
k* = 1, . . . N, where the segment with larger energy has 
smaller slope, thereby implying that the curve is concave. 
The piecewise linear character of microcanonical entropy 
has its origin in the first order poles of the canonical 
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FIG. 2. The energy-temperature relation for N = 3, x ~ 3 f° r 
different L in the (a) microcanonical ensemble and (b) canon- 
ical ensemble. For large values of L, the transition points are 
given by G8x)c = *(**,**- 1). 



partition function, with the slopes being equal to the 
strength of the poles (36j . 

The energy-temperature relation in the microcanon- 
ical ensemble is obtained from the thermodynamic re- 
lation P = || . This gives fix = £(k*,k* - 1), imply- 
ing that /3x takes on N distinct values, corresponding 
to the N transition points, which coincide with those ob- 
tained from the canonical partition function [see Eq. ([5)1], 
In Fig. [2fa), we show the energy-temperature relation 
in the microcanonical ensemble for finite L by consid- 
ering all the terms in g(E). As L is increased we ob- 
tain the limiting step function determined by the rela- 
tion Px = £,{k*,k* — 1). The corresponding data for 
the canonical partition function from Eq. is shown in 
Fig. [2jb). Thus, in the thermodynamic limit (L — > oo), 
we obtain the same limiting curve in both the ensembles. 

We can further quantify the energy-temperature re- 
lation shown in Fig. [2] The energies of the plateaus 
are obtained by solving Eq. ([8j as an equality, and 
are given by E plateau = xN' 1 ^ - k){N - k + 1)L, 
k = 0, . . . , N. A simple physical interpretation can be 
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ascribed to these plateaus, as explained below. Consider 
a scenario when m counterions have condensed onto the 
disc, while the remaining N — m counterions are at the 
boundary R. The energy of this configuration, in the 
limit of large L has contribution from two parts: (1) 
q' 2 L(N — m) 2 corresponding to interaction between a 
disc of charge (N — m)q' and N — m counterions and 
(2) q' 2 L(N — m)(N — m — l)/2 corresponding to interac- 
tion between counterion pairs. The total energy of this 
geometry is thus \N^ 1 (N — m)(N — m + 1). Comparing 
this energy with the energy of a plateau, it is clear that 
the plateau corresponding to a certain value of k corre- 
sponds to a case where k counterions have condensed and 
the remaining ones are at the boundary. 

We stress that we have not used the Legendre trans- 
formation to evaluate the microcanonical entropy from 
the canonical free energy. The entropy obtained by the 
Legendre transformation will always give the concave en- 
velope of the microcanonical entropy [H, |37j ■ Since the 



entropy that we have calculated directly from the den- 
sity of states is concave, we should be able to obtain 
the same by a Legendre transform of the free energy, 
S(E) = (3 E - PF{P), where P is to be eliminated us- 
ing the energy-temperature relation E = A (J3F) . In the 
thermodynamic limit, it is straightforward to do so for 
this model. 

To summarize, we studied counterion condensation 
transition on a two dimensional charged disc in the micro- 
canonical ensemble. In the limit of infinite dilution, we 
obtained an expression for the microcanonical entropy, 
and showed that the entropy-energy curve consisted of 
linear segments with decreasing slope, and hence that the 
entropy is concave with respect to energy. This implies 
the equivalence of the microcanonical and canonical en- 
sembles. In particular, the energy-temperature relation 
and the transition points of CCT obtained from the mi- 
crocanonical entropy are shown to be identical with those 
obtained from the canonical partition function. 
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